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We study ground-state properties of interacting two-component boson gases in a one-dimensional 
harmonic trap by using the exact numerical diagonalization method. Based on numerical solutions 
of many-body Hamiltonians, we calculate the ground-state density distributions in the whole inter- 
action regime for different atomic number ratio, intra- and inter-atomic interactions. For the case 
with equal intra- and inter-atomic interactions, our results clearly display the evolution of density 
distributions from a Bose condensate distribution to a Fermi-like distribution with the increase of 
the repulsive interaction. Particularly, we compare our result in the strong interaction regime to 
the exact result in the infinitely repulsive limit which can be obtained by a generalized Bose-Fermi 
mapping. We also discuss the general case with different intra- and inter-atomic interactions and 
show the rich configurations of the density profiles. 
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I. INTRODUCTION 

The realization of two-component Bose-Einstein con- 
densates (BECs) of trapped alkali atomic clouds [l], [1| 
has stimulated the great interest in the multi-component 
atomic (or spinor) gases. The competition between inter- 
species and intra-species interactions gives rise to a richer 
structure of spinor condensates than their single compo- 
nent counterparts. Many theoretical investigations have 
focused on its stability, phase separation, collective ex- 
citation, dynamics and other macroscopic quantum phe- 
nomena [1, 0) S 01 ■ Moreover resonant control of two- 
body interactions via Feshbach resonances in ultracold 
gases makes it possible to study a two-component con- 
densate with a tunable interspecies interaction. On the 
other hand, the experimental progress in the creation 
and manipulation of ultracold atomic gases in effective 
one-dimensional (ID) waveguides 0, [1,0 offers the op- 
portunity of studying the many-body physics of the ID 
correlated quantum gases even in the strong interaction 
limit. The experimental realization of Tonks-Girardeau 
(TG) gases |lO] in the strongly interacting limit of Bose 
gases has been reported in two recent experiments [1, Q . 

For the ID quantum gases quantum fluctuations 
are greatly enhanced and the effects of strong correlations 
become very important [H, [H, [13, . Although two- 
component quantum gases have been the subject of in- 
tensive theoretical interest 3, 4, 5, 6], most of the studies 
were carried out within the scheme of the mean field the- 
ory which are proven insufficient for strong interactions 
between particles. For the system in the strong interac- 
tion regime, non-perturbation theory has to be exerted. 
For the ID homogeneous system with (5-function interac- 
tion, the two-component Bose gas with spin-independent 
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s-wave scattering (equal interspecies and intra-species in- 
teracting strength) is known to be exactly solvabl e by 
the Bethe-ansatz method in the full physical regime [la] ■ 
The exact results have unveiled that ferromagnetic order 
emerges i n sp inor Bose gases with spin-independent in- 
teraction [3 [13, [3, ■ But for the atomic gases with 
arbitrary interacting strength (with the SU(2) symmetry 
broken) or in an inhomogeneous trap, the system is no 
longer integrable. 

Despite no exact solution for the many-body systems 
with arbitrary interacting strength in a harmonic trap 
available, it is still possible to resort to numerical ex- 
act diagonalization method to obtain the exact solution 
for small size system. Most recently, some sophistical 
theoretical methods, such as the exact diagonalization 
method 2C . l2ll and multi-orbital self-consistent Hartree 
method 23. l23l. [2^ . [25j have been applied to study the 
single-component few-boson systems in a harmonic trap. 
Due to properly accounting the many-body effects, such 
methods work well even in strongly interacting regime 
and yield consistent results with the exact results ob- 
tained by the Bethe-ansatz method [2^, [l^l . In this work, 
we shall study two-component boson systems with repul- 
sive interactions and focus on few-atom systems which 
can be treated exactly by exact diagonalization method. 
One of the purpose of the present work is to address the 
role of the correlation in the trapped interacting two- 
component system under a well controllable way, in par- 
ticular in the strongly interacting limit where the approx- 
imate mean-field methods are generally not able to give 
correct results. Our study can display the transition from 
a weakly interacting condensate to the strongly interact- 
ing limit in which the strongly repulsive interactions be- 
tween atoms plays a role of effective Pauli principle. In 
this limit, we compare our result with the exact result 
in TG-gas limit via a generalized Bose-Fermi mapping 
which has been applied to the Bose-Fermi mixtures by 
Girardeau et. al very recently We also address the 
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case with different inter-species and intra-species inter- 
actions. 

Our paper is organized as follows. Section II introduces 
the model and gives a brief introduction to the compu- 
tational method. The subsequent section is devoted to 
our results for the ground-state density distribution. A 
summary is given in the last section. 



II. THE MICROSCOPIC MODEL AND THE 
METHOD 

Let us consider the two components of Bose gases con- 
fined in a strongly anisotropic harmonic trap, where the 
transverse confinement frequency ujj_ is strong enough so 
that the dynamics of Bose gases are effectively described 
by a ID model due to the radial degrees of freedom are 
frozen. We start with the microscopic Hamiltonian of 
two species of Bosons 



harmonic oscillator 



H = 



52 



+gi2 / dx^\{x)^\{x)^2{x)^i{x) 



+ Va,{x) 



^c.{x) 



(1) 



with Va{x) 



(a = 1, 2) being the trap poten- 



tial felt by the a-component bosons with mass nia and 
longitudinal trap frequency lOq, . Here vPq, {x) denotes the 
field operator of the a-component bosons for the longitu- 
dinal degree of freedom. In the following we will focus on 
the two-component Bose gases composed of two hyperfine 
state of the same atoms such that mi = 7712 = m and it is 
assumed that both of them feel the same trap frequency, 
i.e., ui — LJ2 = ^- The effective intra- and inter-atomic 
interaction strengths are denoted as ga = ~2?i^/(ma"^) 
(a = 1,2) and 512 — —2h^/{ma\^), respectively. They 
are related to the three dimensional s-wave scattering 
length between two atoms (aj and being the intra- 
atomic scattering len gth and al^ being the inter-atomic 
scattering length) by [111 . fl2| 



IC(1/2)| 
V2 



a = 1,2,12 



with lj_ — y^h/muj^ the transverse characteristic oscilla- 
tor length. 

Expanding the operator '^a (x) in the basis of harmonic 
oscillator functions (orbital) 4>i (x) with i the orbital index 



i=0 



we get the many-body Hamiltonian in terms of the cre- 
ation and destruction operators bl^ and bia of the axial 



2 ^ U,jMOiJ)]j}kabla 

i,j,k,l 



+ U12 ^ Ii,],k,lb\J}]2^k2bll 

i,j,k,l 



(2) 



with fii = i + i. Here we have rescaled the energy and 

length in units of huj and Iq = y^h/muj. The parameters 
Ua {a — 1, 2, 12) are dimensionless interaction constants 
and li-jki = lo J^^dx(l)i{x)(j)j{x)(f>kix)(l)i{x) are dimen- 
sionless integrals. The ground state of the Hamiltonian 
can be obtained by diagonalizing it in the Hilbert space 
spanned by the eigenstates of the harmonic oscillator. 
Although the exact diagonalization methods can in prin- 
ciple solve the ground-state problem of the many-body 
Hamiltonian exactly, one has to truncate the set of ba- 
sis functions to the lowest L orbitals for reasons of nu- 
merical feasibility. The dimension of the space should be 
^ ^N2+L-i if -^1 1-componcnt atoms and N2 2- 
component atoms are populated on L orbitals. In princi- 
ple the exact ground state of the many-body Hamiltonian 
can be obtained even in the strongly interacting regime as 
long as L is large enough. Our goal is to investigate the 
ground states of the interacting system for all relevant 
interaction strengths in a numerically exact way which 
however restricts our study to small particle system. We 
note that either the effort to extend the exact diagonal- 
ization method to solve the larger system or to get a 
well controllable numerical result in the strongly inter- 
acting limit is a highly challenging and time-consuming 
task as what have been displayed in the works for the 
single-component model [20, [23, [23] • 



III. GROUND STATE PROPERTIES 

In the following we will investigate both the density 
distribution of the many-body ground state of each com- 
ponent and the total distribution, which is given by 



P{x) = ^"(^) 



a=l 



with 



Po,{x) = (GSI^-Kx) *a(a;)|GS) 

- 5]0*(a;)<^,(a;)(GS|SL&,a|GS). (3) 

The ground state |GS) can be obtained by numerically di- 
agonalizing the Hamiltonian in the reduced Hilbert space 
which is typically composed of 10^ basis states in our cal- 
culation. To guarantee the numerical exactness, in the 
following we will focus on the system of four bosons with 
the maximum energy cutoff of SOtuv corresponding to the 
number of orbitals L=30. 
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FIG. 1: Density distribution of the ground state for A'^i — 
N2 = 2. (a) [/i = U2 = U12 = 0.1; (b) Ui = U2 = U12 = 3.0; 
(c) Ui=U2 = U12 = 7.0; (d) Ui = U2 = U12 = 10.0. 



FIG. 2: Density distribution of the ground state for A''i = 1 
and iV2 = 3. (a) Ui ^ U2 = U12 = 1.0; (b) Ui ^ U2 ^ U12 = 
3.0; (c) Ui = U2 = U12 = 5.0; (d) f/i = C/2 = U12 = 10.0. 



We first consider the case with the intra- and inter- 
atomic interactions being same, i.e., Ui = U2 = U12. In 
Fig. 1, we display the density distributions of two com- 
ponents of Bose gas in a harmonic trap when both com- 
ponents have same particle numbers. As shown in the 
figures each component of the Bose gas has the same dis- 
tribution in the whole regime of interaction. In the weak 
interaction regime, the Bose gas exhibits a Gauss-like 
distribution, which implies that the many-boson system 
can be approximately described by a single wave func- 
tion, i.e., all the bosons condense into the same single 
state. In the scheme of the mean-field theory, the sin- 
gle wave function is obtained by solving Gross-Pitaevskii 
equations. In the limit of [/„ = (a = 1, 2, 12) all atoms 
occupy the lowest orbital and the ground state takes the 
form of 

1=1 3 = 1 

which can be expressed as |GS) = (&oi)^H&o2)^^ |0) ™ 
our Hilbert space. With the increase of the atomic inter- 
action, the density distribution becomes broader and flat- 
ter. When it reaches the strongly interacting regime, four 
peaks appear at Ua — 10 as shown in the Fig. Id. The 
emergence of density oscillations is a signal of "fermion- 
ization", i.e., the density profile shows fermion-like dis- 
tribution [l^iUll. In this case, two bosons would avoid to 
occupying the same orbital and thus the bosons occupy 
the orbitals one by one similar to the fermions. Therefore 



the mean-field theory based on the single wave function 
approximation is not a suitable theory in the strongly in- 
teracting limit. To capture the character of density oscil- 
lations, more sophistical multi-orbital mean-field theory 
is required [12, [ll, [H, [11] . 

For the imbalanced case with A^i 7^ N2, the total 
density distribution displays the similar crossover be- 
havior as the case with A^i = iV2. As shown in Fig. 2, 
the halfwidth of the total density distributions becomes 
larger and larger with the increase of interaction strength, 
followed by the emergence of an obvious spatial oscilla- 
tion structure at Ua = 10. What special here is the den- 
sity distribution for each component is no longer same 
in the whole interacting regime. In the weakly interact- 
ing regime, each component has similar density profile. 
While in the middle regime, the component of the Bose 
mixture behaves in different way. At Ua = 3, two obvious 
peaks are observed for the component with less atoms 
whereas the total density shows no oscillation. With 
further increasing the interaction strength to Ua = 5, 
the total density distribution begins to oscillate with two 
peaks emerging. When reaching the strongly interacting 
regime with Ua — 10, four obvious peaks are observed 
for both components of the Bose mixture. In this limit 
each component has similar density profile again: the 
density distribution of the component with Na atoms 
is Na/{Ni + N2) of the total density distribution with 
Ni + N2 peaks. As we shall discuss latter, this is con- 
sistent with the result in the limit of infinite interaction, 
where both the total distribution of two species and the 



4 




FIG. 3: Density distribution of the ground state for A''i — 
N2=2 and Ui = U2 = 0.1. (a) U12 = 0.1; (b) U12 = 3.0; (c) 
U12 = 5.0; (d) U12 = 10.0. 

respective distribution of them behave hke noninteract- 
ing spinless fermions. 

In the hmit of = 00 , we can app ly the generahzed 
Bose-Fernii mapping method [28l . [23 | to construct the 
ground state wavefunction of the Bose mixture. Accord- 
ing to Girardeau's argument [lo| . in the infinite interact- 
ing hmit the interactions can be replaced by constraints 
that the many-body wave function vanishes at all colli- 
sion points [IoIjIIIjIIOI- It follows that the Girardeau- type 
ground state wavefunction has the following form 

^/'(a;]^ , . . . , Xj^_^ ;Xi, . . . ^ Xj^^ ) 

— -^{^1 I ■ ■ • I ^ATj T^lT ■ ■ ■ , x^^ ) X 

ipF{xl,...,xlf^;xl,...,xjjJ (4) 

with the antisymmetric function A(xi, X2, ■ ■ ■ , xjy) — 
ni<ij<ArSgn(a:, - Xj) and ippixl, . . . , x^^^- xj, . . . , x%J 
is a Slater determinant of N — Ni + N2 orthonormal or- 
bitals (poix), • • • , (t>N-i{x) occupied by both components 
of bosons. Explicitly, we have 

...,xn) = (^!)-^det[0,(xO]^r°j;;;;;iJ^-' 

with xi , . . . , xn denoting x\, . . . , xl^_^ ]x\, . . . , . Now 
it is straightforward to get 

Pa[x) = ^pf{x) (5) 

with Pf{x) — X]il;0^ I '/'i (2^)1^ : i-6., the densities of both 
components are proportional to the density of a harmon- 
ically trapped TG gas of N bosons which is identical to 



FIG. 4: Density distribution of the ground state for A'^i = 1, 
TVs = 3 and ^ U2 ^ 0.1. (a) U12 = 0.1; (b) U12 = 1.0; (c) 
U12 = 2.0; (d) U12 = 10.0. 



the density of N free fermions. In Fig. 2d, we compare 
our numerical result with the exact result in the TG limit. 
It is clear that the numerical result complies with the ex- 
act results very well which implies that the system has 
entered the strongly interacting regime at this parame- 
ter point. Our numerical results provide a clear evidence 
that the state obtained by Girardeau-type construction 
|28t remains to be a good approximation for large but 
finite repulsion. 

Now we discuss the general case with different intra- 
and inter-atomic interactions. First we calculate the den- 
sity distributions for the equal-mixing Bose mixture with 
equal weak intra-component repulsions. For convenience, 
we take Ux — U2 — 0.1 and Uu ^ Ui,U2. Density dis- 
tributions for different inter-component interactions are 
shown in Fig. 3. In all the interaction regime, we have 
Pi{x) — p2{x) — \p{x). With the increase of C/12, the 
density distribution becomes broader. For U12 — 5, two 
peaks occur. With further increase of inter-component 
interaction, only two peaks in the density distribution 
are observed even for U12 — 10. Due to the intra-atomic 
interaction is weak, each component of atoms may form 
a condensate. Similar results are found for the system 
with weak intra-component repulsions even ii Ui 7^ U2- 
The difference from the plot of the case of Ui = U2 is 
that the component with stronger intra-atomic interac- 
tion distributes in broader regime than that with weaker 
interaction. 

For the case of Ni ^ N2 with weak intra-atomic inter- 
action Ui = U2 = 0.1, density distributions for different 
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inter-component interactions are displayed in the Fig. 4. 
Similar to the equal-mixing case, the total density dis- 
tribution becomes flatter with the increase of C/12. Al- 
though two obvious peaks occur for the component with 
less atoms when U12 — 2.0, no obvious peak is detected 
in the total density profiles even for U12 = 10.0. When 
the inter-atomic interaction is strong, the density for the 
different component has different spacial distribution and 
gets demixed in order to lower the inter-species interac- 
tion energy. As shown in the figures, the component with 
more atoms stays in the middle of the trap and the com- 
ponent with less atoms stays in the regime away from the 
middle of the trap. 

Finally, we study the change of density distributions 
versus different intcr-componcnt interactions for the Bose 
mixture with strong intra-atomic interactions. For con- 
venience, we consider the imbalanced case with iVi = 3 
and 7V2 = 1 and we take equal strong intra-component 
repulsions Ui = U2 = 10.0. For the vanishing inter- 
component interactions, the atoms of two species al- 
ready lie in the strongly interacting regime and the den- 
sity distribution for each component gets the Fermi-like 
feature. The distribution of pi{x) displays three peaks 
and the distribution of P2(x) displays one peak. As 
the inter-atomic interaction increases, the distribution of 
each component changes continuously as shown in Fig. 
5. When U12 — 3.0, the central peak of the first compo- 
nent is suppressed and the density profiles get broader. 
At U12 — 8.0, four peaks occur for the first component 
whereas the density for the second component get more 
flatter accompanying with two peaks discernable. With 
further increase to U12 = 10.0, both the total density dis- 
tribution and the distribution of each component show 
Ni + N2 peaks and the latter takes on the same distri- 
bution normalizing to their own particle numbers N^. 



IV. SUMMARY 

In summary, we have investigated the ground-state 
properties of two components of Bose gases composed 
of two hypcrfinc states of the atoms. With the nu- 
merical diagonalization method the density distributions 
are obtained, which display obviously different properties 
for different parameters in the system. For the system 
with equal intra-component and inter-component inter- 
actions, the total density distribution shows similar prop- 
erties to the single-component Bose gases with Ni + N2 
atoms. The evolution from weakly interacting conden- 
sate to strongly interacting Tonks gas is shown with the 
increasing atomic interaction. In the weakly interact- 
ing regime the Bose mixture displays Gauss-like distri- 
bution whereas in the strongly interacting regime the 
density distribution gets Fermi-like features. Our nu- 



merical result is consistent with the exact result in the 
infinite limit and thus provides numerical evidence that 
the Bose mixture can be approximately described by the 
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FIG. 5: Density distribution of the ground state for Ni = 3, 
Ni = 1 and (7i = f/2 = 10.0. (a) U12 = 0.0; (b) U12 = 3.0; (c) 
U12 = 8.0; (d) U12 = 10.0. 



TG gas even for the system with strong but finite repul- 
sion. We also discuss the general case with different intra- 
component and inter-component interactions and display 
the evolution of density distributions from weak to strong 
inter-component interaction regime for both cases with 
weak and strong intra-component interactions. Despite 
the size limitation of the exact diagonalization method 
used in the present work, our results are no doubt very 
meaningful because it provides us important insights into 
the quantum many-body physics beyond various approx- 
imate effective-field approaches. 

Note added: After finishing the present work [3l| . 
we noticed a related work by ZoUner et. al [32*1, in 
which similar issue was discussed by the multi-orbital 
self-consistent Hartree method. 
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